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ABSTRACT 

Weak gravitational lensing is a promising probe of dark matter and dark energy requiring 
accurate measurement of the shapes of faint, distant galaxies. Such measures are hindered by 
the finite resolution and pixel scale of typical cameras. On the other hand, as imaging telescopes 
are practically limited to a fixed number of pixels and operational life-span, so the survey area 
increases with pixel size. We investigate the optimum choice of pixel scale in this trade-off for a 
space-based mission, using the full engineering model and survey strategy of the proposed Super- 
Nova/ Acceleration Probe as an example. Our methodology is to simulate realistic astronomical 
images of known shear and to evaluate the surface density of sources where the shear is accurately 
recovered using the Rhodes, Refregier & Groth algorithm in the context of the derived dark mat- 
ter power spectrum. In addition to considering single exposures, we also examine the benefits of 
sub-pixel dithering. Although some of our results depend upon the adopted shape measurement 
method, the relative trends, particularly those involving the surface density of resolved galaxies, 
are robust. Our approach provides a practical counterpart to studies which consider the effects 
of pixelation from analytic principles, which necessarily assume an idealized shape measurement 
method. We find that the statistical error on the mass power spectrum is minimized with a 
pixel scale equal to 75-80% of the FWHM of the point-spread function, and that dithering is 
marginally beneficial at larger pixel scales. 

Subject headings: gravitational lensing — surveys — instrumentation: detectors 



1. Introduction 

The gravitational field arising from a massive 
foreground structure deflects and distorts the light 
from distant objects according to the process of 
gravitational lensing (c/ Blandford et al. 1991; 
Narayan & Bartelmann 1996). Rapid progress has 
been made in the accurate measurement of the 
weakly distorted shapes of background galaxies in 
order to determine the distribution of dark matter 
in space and time (for reviews see Mellier 1999; 
Bartelmann & Schneider 2001; Refregier 2003b). 
As the growth of the dark matter power spectrum 
is also sensitive to both the gravitational attrac- 
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tion of dark matter and the repulsive effect of dark 
energy, weak gravitational lensing is emerging as a 
very promising cosmological probe (Dark Energy 
Task Force, Albrecht et al. 2006). 

Several proposed space-based missions (SNAP, 
Aldering et al. 2004; DUNE, Refregier et al. 2006) 
plan to use weak lensing as a major probe. A 
key design element is the detector pixel scale, # cc d 
(arcseconds), which must be optimized in terms 
of a trade-off between the fidelity with which the 
true weak lensing signal is recovered and the ben- 
efits of a large survey area during a fixed mission 
lifetime. In this paper, we investigate the prac- 
tical consequences of changing the detector pixel 
scale, using simulations to explore empirically the 
effect on measurements of the cosmological matter 
power spectrum. 

Extracting the weak lensing signal requires ac- 
curate measurements of the shapes of faint galax- 
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ies, which are inevitably degraded by the convo- 
lutions arising from a finite point spread function 
(PSF) and detector pixelation. The PSF of an op- 
tical space telescope will be dominated by diffrac- 
tion arising from the system pupils, mirror struts, 
and other physical structures which perturb inci- 
dent light. The PSF blurs the shapes of back- 
ground galaxies, which need to be deconvolved or 
at least corrected. The shape of the PSF, which 
typically varies as a function of detector position 
and time can normally be obtained by examining 
point sources or directly by ray-tracing software. 

Pixelation similarly degrades faint image shapes. 
In the case of a CCD detector, a photon incident 
at a photosite, or pixel, photoelectrically liber- 
ates an electron, which is then bound near the 
CCD surface to be measured later. The photosite 
has finite cross-sectional area and counts photons 
that fall anywhere inside its perimeter, averaging 
over all sub-pixel scale image features. Liberated 
electrons may also diffuse to the surface of adja- 
cent photosites, an effect termed charge diffusion. 
Averaging and diffusion will distort the measured 
shape of features nearly the size of a pixel or 
smaller, biasing the measured weak lensing signal. 

While CCD photosites effectively average over 
small scale image information, it is nonetheless 
possible to recover some of the high frequency in- 
formation by dithering. Dithering involves tak- 
ing multiple exposures of the same stationary ob- 
jects so that their sub-pixel positions in the CCD 
are different each time. This can be done by 
slightly translating the camera between exposures. 
Individual dithered exposures are combined be- 
fore shape analysis using the Drizzle algorithm 
(Fruchter & Hook 2002) or the Fourier techniques 
of Lauer (1999a, b); or during shape analysis by 
somehow averaging the measurements. 

The Drizzle algorithm averages dithers on a 
pixel grid that can be finer than $ C cdi 

resulting 

in a higher resolution, resampled image. Lauer's 
technique accomplishes the same thing without re- 
pixclation by combining Fourier transforms from 
the individual dithered exposures. We investigate 
the extent to which Drizzling can recover weak 
lensing signals from undersampled data. 

Our results are part of a wider investigation of 
the accuracy of weak lensing measurements and 
the effects of various instrumental and algorithmic 
parameters. Our image simulation software is em- 



ployed by the Shear TEsting Program 1 (STEP). 
The STEP project comprises independent "blind" 
analyses of simulated data of various types by in- 
ternational groups with a goal of verifying the lim- 
itations of extant weak lensing measurement al- 
gorithms. Simulations of both ground (Heymans 
et al. 2006; Massey et al. 2007b) and space-based 
(in prep.) data have been undertaken, and in par- 
ticular space STEP will also explore the effects of 
pixelation on shear recovery. Bernstein (2002) and 
Bernstein et al. (in preparation) perform a simi- 
lar analysis, but starting from the opposite end of 
analytic first principles, studying the irreversible 
loss of information during pixelation, even with a 
perfect shape measurement method. In a further 
study, Massey et al. (2004b) analyze the effect of 
varying the exposure time on weak shear recov- 
ery from a SNAP-like mission. Refregier et al. 
(2004) study semi-analytically the trade-off be- 
tween a wide and deep SNAP survey strategy by 
looking at the expected errors on cosmological pa- 
rameters from weak lensing data. 

This paper is organized as follows. In §2 we de- 
scribe a useful mathematical model of pixelation. 
§3 introduces our adopted PSF model (§3.1) and 
simulation software (§3.2) and discusses how we 
vary the pixel scale (§3.3), our simulated dither 
and Drizzle strategy (§3.4), and our shear mea- 
surement algorithm (§3.5). Results are presented 
in §4. Properties of the simulated images are 
shown in §4.1, and then the recovered weak lensing 
signals are analyzed in terms of the surface den- 
sity of useable galaxies (§4.2), the bias determi- 
nation (§4.3), the standard error of shear estima- 
tion (§4.4), and the predicted error on the matter 
power spectrum (§4.5). We conclude in §5. 

2. A Pixelation Model 

Here we introduce the formalism necessary for 
our analyses. We define the true time- independent 
image I tT (6) as the monochromatic photon flux at 
angular location 6 = (9 X , 6 y ) in the sky, as imaged 
with a perfect instrument. A real telescope con- 
volves I tr (6) with its "diffraction-pattern" PSF, 
P(#), which is normalized to unit integral. In our 
simulations, we take the PSF to be spatially inde- 
pendent across the focal plane. 



1 See http://www.physics.ubc.ca/~heymans/step.html. 
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The PSF-convolved image is then pixelated. In- 
tegration within a pixel is equivalent to two dis- 
tinct operations (Lauer 1999a, b). First, the PSF- 
convolved image is again convolved with a pixel 
response function (PRF) R(0), which has charac- 
teristic size 0r, producing 



I 1 (6) =/ tr (0)*P(0)*R(6/), 



(1) 



where * denotes convolution. Second, I 1 is sam- 
pled at regular intervals of size 9 Sl such that the 
observed ij-th pixel value is 



(2) 



where i G {1,2,..., N x }, j e {1,2,..., N y }, and 
the observed image has N x x N y pixels. 

It is important to emphasize that pixelation is 
more complicated than a re-sampling of the PSF- 
convolved image. Photosites integrate photons in- 
cident anywhere within their small yet finite bor- 
ders, allowing only for a spatially averaged, or 
"binned," flux measurement. The PRF convolu- 
tion of Equation [1] fully describes this averaging. 
CCDs become ideal samplers only in the limit of 
truly infinitesimal photosites, in which case the 
PRF is a 5-function and the convolution is simply 
equivalent to the PSF-convolved image. 

The ideal CCD with finite photosite size has a 
single PRF across the entire array, equal to a unit- 
normalized square tophat with side 6*r = 9 cc d- 



Rsq(#) — 



Cd if \U\0y\ <0 ccd /2, 







otherwise. 



(3) 



Real PRF's deviate from this idealized case. For 
example, imperfect quantum efficiency, which re- 
sults in only a fraction of incident photons being 
counted, is equivalent to an R normalized to less 
than unity. Similarly, variable quantum efficiency 
across the array is equivalent to a spatially- varying 
R; and charge diffusion can be regarded as a PRF 
with 6> R > 9 ccd . 

We can apply this formalism to consider the 
performance of single exposures and dithered im- 
ages ("dithers") of various kinds. A single expo- 
sure is an image Iij with sample spacing 9 S — 9 cc d- 
A dither refers to the case where the camera is 
translated by, e.g. non-integral pixels (d x ,d y ) de- 
livering an observed image Ii_d x ,j-dy, also with 
9 S = #ccd- Ideal interlacing is a particularly im- 
portant case consisting of N x N dithers with 



d x = k/N and d y — £/N, where N is a posi- 
tive integer, and k,£ G {0, 1, . . . , iV — 1}. Driz- 
zle can dcintcrlace such dithers, say from half- 
$ccd shifts (N = 2), by applying them to a pixel 
grid that is twice as fine, yielding an observed im- 
age with 9 S = #ccd/2. Drizzling in general is an 
effective "re-pixelation," with its own PRF that 
is free to be chosen by the user. We emphasize 
that the observed image in all cases is convolved 
with the CCD pixel kernel R with size and 
9r 7^ #ccd 7^ 9 S in general. That is, even in the 
limit of infinite, perfectly deinterlaced dithers with 
a 5-function re-pixelation kernel, the resultant im- 
age is still convolved with R. 

We model the PRF as a square top-hat response 
convolved with an additional charge diffusion ker- 
nel, multiplied by a fraction e q e (0, 1] represent- 
ing quantum efficiency. The charge diffusion ker- 
nel is taken to be a Gaussian, 



Rch(#) 



1 



cxp 



(4) 



where a c h is the RMS extent of the diffusion. Our 
PRF is then 



R(0) 



e q R sq (0)*R ch (0). 



(5) 



The final size of the PRF, 9r, is larger than both 
# cc d and CT c h due to the convolution. 

Following Bernstein (2002) and assuming P and 
R are both space- and time- invariant, we use the 
associativity of convolutions to define the effective 
point-spread function ePSF, equal to P * R. Now 
I tT (0) is convolved only once, with the ePSF, and 
then sampled as in Equation 2. Throughout the 
paper we distinguish between the PSF, P, which 
only includes diffraction; the chPSF, equal to P * 
R c h; and the ePSF, equal to P * R c h * R sq - 

3. Method 

We have developed a pipeline to simulate all rel- 
evant steps from the acquisition of imaging data, 
through its reduction, to the measurement of the 
weak lensing (or "shear") signal. We first manu- 
facture realistic images containing galaxies with a 
known shear signal, with various pixel scales, but 
keeping all other parameters fixed. We then de- 
tect galaxies in the noisy images and measure their 
shapes. We correct their shapes for ePSF effects, 
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which are measured from separate, simulated im- 
ages of dense stellar fields at each pixel scale. We 
finally compare the output, recovered shear mea- 
surement to the known input shear signals. In a 
parallel pipeline, we also create sets of four shal- 
lower but dithered images at slightly larger pixel 
scale, which we stack using Drizzle to improve the 
pixel sampling. The interlacing of these images 
is optimal, thus providing an optimistic study of 
how much resolution could in principle be recov- 
ered from a hardware design that slightly under- 
samples the PSF. 

In the following subsections we describe the key 
ingredients and processes in the pipeline. 

3.1. Adopted PSF 

The assumption of a constant ePSF in space or 
time isolates the problem of shear measurement 
(in which we are primarily interested) without dis- 
tractions of PSF interpolation (which is a separa- 
ble problem that is being widely discussed else- 
where) . 

Our ePSF model is based on an early en- 
gineering design, called TMA63, of the Super- 
Nova/ Acceleration Probe 2 space telescope (SNAP, 
Aldering et al. 2004), and accounts for diffraction 
by the telescope plus charge diffusion in the CCD. 
The simulated diffraction pattern comes from a 
raytraced model of the SNAP //ll optical sys- 
tem. It simulates light from an 820nm wavelength 
point source incident on a 2m primary mirror, a 
secondary supported in front of the primary by 
three struts, a folding flat, and a tertiary. The 
point source is taken to be off-axis such that its 
image appears on the focal plane 198mm radi- 
ally away from the optical axis, where the typical 
SNAP CCD lies. SNAP plans to employ Lawrence 
Berkeley National Laboratory's new high resis- 
tivity CCDs, in which charge diffusion further 
spreads point source light (Holland et al. 2003; 
Stover et al. 2000; Groom et al. 2000). We model 
charge diffusion as an additional Gaussian convo- 
lution with standard deviation cr c h = 4^im, the 
expected level for these CCDs, which have 10/im 
wide photosites. The net chPSF has FWHM of 
0.12"; pixclation at the baseline SNAP CCD pixel 
scale of 0.10" gives a final ePSF FWHM of 0.14". 



2 See http://snap.lbl.gov/. 



3.2. Simulated Images and Input Shear 

The image simulation suite of Massey et al. 
(2004a), Simage, is our main tool for creating ar- 
tificial astronomical images and applying an ar- 
bitrary weak shear signal. Simage uses Shape- 
lets 3 , a parametrization of galaxy morphologies as 
a weighted sum of a complete, orthonormal set of 
basis functions (Refregier 2003a). 

Realistic morphologies are generated by empir- 
ically matching the measured properties of ac- 
tual galaxies in the Hubble Deep Fields (HDF, 
Williams et al. 1996, 1998). The HDF galaxy 
positions, orientations, morphologies (plus sizes 
and magnitudes) are randomized when generat- 
ing new images. In this way realistic magnitude- 
morphology trends are produced, though no spa- 
tial clustering is imposed and no redshift informa- 
tion is encoded . Because the HDF source cat- 
alog on which the simulations are based are al- 
ready pre-convolved with the HDF PSF, they have 
slightly larger intrinsic sizes than the true galaxy 
population. However, this does not affect the pro- 
cess of shear addition and measurement. 

The analytic shapelet models of galaxies are 
then sheared as described by Massey & Refregier 
(2005), including terms up to fourth order in the 
shear 7. Our input shears 7 = (71, 72) range from 
—0.05 to +0.05 in steps of 0.025 in each compo- 
nent, while the other component is fixed at 0, viz: 

71 e {-0.05, -0.025, 0, 0.025, 0.05} and 72 = 0; (6) 
71 = and 72 e {-0.05, -0.025, 0.025, 0.05}. (7) 

A total of nine different input shears per pixel 
scale is thus produced. We make one image per 
input shear, uniformly distorting all galaxies in 
the field. The sheared galaxies are convolved with 
the PSF and transformed to real image space, in- 
tegrating the basis functions analytically within 
each pixel. This integration is mathematically 
identical to the convolution with the R sq described 
in §2. Finally, photon noise and a realistic space 
background signal are added. 



3 See http://www.astro. caltech.edu/~rjm/shapelets/ for a 
Shapelet analysis package. 
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3.3. Pixel Scale 

The key question motivating this study is what 
pixel scale most benefits weak lensing analyses. To 
address this we simulate and analyze images with 
CCD pixel scales perturbed about a baseline value 
chosen to be 9 — 0.10". In practice, pixel scales 
can be adjusted in one of two ways and this defines 
two categories of simulations. 

First, we simulate changing the CCD pixel scale 
$ccd by changing only the physical size of the pho- 
tosites, labeling this dataset "Sim 1." Second, we 
simulate changing cc( j by adjusting only the focal 
length of the telescope, which changes the plate 
scale, calling this dataset "Sim 2." While these 
each linearly perturb 9 cc d, their detailed effects on 
the final image are distinct and therefore each mer- 
its analysis. 

The Sim 1 dataset consists of different simu- 
lated images with detector pixel scales between 
0.04" and 0.16" at 0.01" intervals. This ranges 
from the resolution of the processed, end-product 
HDF images, to just larger than the baseline ch- 
PSF size. We adopted a baseline photosite size of 
10/im, and the different pixel scales are achieved 
by varying this between 4/j,m and 16^im in 1/xm 
steps. In practice, such changes would need to 
be made at the CCD manufacturing stage. Both 
the diffraction pattern and the baseline plate scale 
of 0.01"/^m are unaffected. Thus, assuming the 
electron diffusion length is constantly 4/im inde- 
pendent of the photosite size, the different photo- 
site sizes cause no net effect in the apparent an- 
gular size of the charge diffusion. The final shape 
and angular size of the chPSF, that is before pixc- 
lation, are left unchanged at all # cc d- This scheme 
for pixel scale adjustment represents a relatively 
academic exercise into the effect of pixelation of a 
fixed PSF. 

The Sim 2 dataset emulates the effect of chang- 
ing the focal length while keeping the physical 
size of the CCD photosite fixed, thereby chang- 
ing the plate scale and 9 cc d from 0.04", to 0.16". 
At longer focal lengths, the Airy disc from diffrac- 
tion, which is constant in angular size, grows lin- 
early in microns at the focal plane while the pixel 
scale decreases in arcseconds and the charge diffu- 
sion remains at 4/^m. At these small pixel scales, 
the diffraction pattern dominates the chPSF as 
the first ring is readily apparent. At shorter fo- 



cal lengths, the Airy disc shrinks while the pixel 
scale increases. In this case the charge diffu- 
sion dominates, smoothing out the diffraction fea- 
tures and broadening the chPSF. This is perhaps 
a more practical engineering solution to adjusting 
the pixel scale. 

In both cases, the 9 cc d = 0.04" images contain 
4096 x 4096 pixels. These dimensions change lin- 
early with 9 cc d such that the subtended solid an- 
gle per image is constant. This does not neces- 
sarily simulate how an instrument would change 
its CCD array size with 9 cc d in practice. We add 
about 6000 galaxies and no stars to each sheared 
image, extending their distribution well below the 
intended detection threshold. Adding noise to fix 
the survey depth to m ~ 27.7 with 2000s expo- 
sures, we reproduce the number of galaxies useful 
for weak shear estimation, n g w 100 per arcmin 2 , 
found at the baseline pixel scale in earlier studies 
(Massey et al. 2004b). For each pixel scale, we also 
make one additional image containing only stars, 
simulating a weak lensing survey that periodically 
points at stellar fields to characterize the ch- and 
cPSF. 

Figure 1 illustrates the effect of perturbed CCD 
pixel scale on the chPSF profile in the two simu- 
lation sets. The net result is that the underlying 
chPSF in all Sim 1 images is constant, while in 
Sim 2 the chPSF shape and size change with # C cd- 
Pixelation further distorts the chPSF to produce 
the ePSF, which we measure and plot in Figure 
2. Also plotted in Figure 2 are ePSF sizes after 
Drizzling, which we discuss presently. 

3.4. Simulating and Combining Dithers 

In addition to the single exposures introduced 
above, we also create two additional datasets, la- 
beled "Sim 1 Drizzled" and "Sim 2 Drizzled" , 
which simulate dithers that we combine to explore 
the extent to which resolution can be recovered, as 
described in §1 and §2. Sim 1 Drizzled emulates 
Sim 1 in how 9 cc d varies by changing the size of 
CCD photosites, and Sim 2 Drizzled emulates how 
Sim 2 changes the pixel scale by changing the plate 
scale. 

In the Drizzled datasets we implement the sim- 
plest, ideally interlaced 2x2 dither pattern, where 
four exposures are taken, each shifted by exactly 
half a pixel in orthogonal directions. These dithers 
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-0.2 -0.1 0.0 0.1 0.2 

arcsec 



Fig. 1. — Cross sections through the peak of 
the simulated point-spread function on linear (top 
panel) and log (bottom panel) scales, with compo- 
nents from diffraction and charge diffusion, then 
finely pixelated at 0.01"; this is what we term 
the chPSF. The baseline chPSF (black line) comes 
from raytracing simulations of the current SNAP 
design, convolved with a 2D Gaussian to approx- 
imate charge diffusion at the CCD. In this study 
we perturb the baseline SNAP detector pixel scale 
of 0.10" in two realistic ways, in parallel, which 
have distinct effects on the chPSF shape. First, 
we change only the CCD photosite size in microns, 
which keeps the angular size of the charge diffu- 
sion kernel constant and therefore leaves the pre- 
pixelated PSF shape unchanged. Second, we ad- 
just only the focal length, which changes the plate 
scale and thus the angular size of the charge diffu- 
sion. In this case, a smaller plate scale reduces the 
effect of charge diffusion over the diffraction pat- 
tern in the final chPSF (light gray line), while a 
shorter focal length causes the diffraction pattern 
to dominate (light gray line). 
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Fig. 2. — Measured full- width- at-half-maximum 
values of the ch- and cPSFs. The chPSF in the Sim 
1 (Drizzled) is independent of the CCD pixel scale, 
while the chPSF in the Sim 2 (Drizzled) increases 
in angular size with increased 9 cc d- The plotted 
chPSF FWHM values (dash lines) are measured 
from highly oversampled images using Source Ex- 
tractor. Shapclct coefficients are measured from 
the oversampled images, and later re-pixelated in 
the simulated astronomical images at the plotted 
CCD pixel scales. The ePSF FWHM's shown are 
averaged from Source Extracted images contain- 
ing only stars. The Drizzled images are resampled 
on pixel grids that are twice as fine as # cc d; we 
measure stars on the Drizzled images, though we 
plot their CCD pixel scale values instead to allow 
comparison with their true underlying chPSFs. As 
this plot shows, pixelation increases the apparent 
size of the point-spread function due to the convo- 
lution. Drizzling introduces an additional pixela- 
tion, which further increases the ePSF size. These 
effects must be considered in precise shape analy- 
ses such as weak lensing requires. 
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resample the ePSF-convolved image from Equa- 
tion 1 at shifted intevals (d x ,d y ), where 



Dither 1: 


(d x , 


dy) 


= (0,0) 


(8) 


Dither 2: 


{d x , 


dy) 


= (0.5,0) 


(9) 


Dither 3: 


(d x , 


dy) 


= (0.5,0.5) 


(10) 


Dither 4: 


(d x , 


dy) 


= (0,0.5) 


(11) 



We then dcinterlace them using Drizzle, placing 
the four dithers on a pixel grid that is twice as 
fine, such that 9 S = # cc( j/2, and setting Drizzle pa- 
rameter pixfrac=0.5. We map each dither pixel 
to one Drizzled pixel with no overlap, such that 
the noise between adjacent pixels is uncorrelated. 
Note that this is an idealized case, impossible to 
achieve in practice, and therefore represents an op- 
timistic case study. 

All dithers and Drizzled images subtend the 
same solid angle as in the other two studies. Each 
dither is a 500s exposure, resulting in a 2000s ef- 
fective exposure after Drizzling. The CCD pixel 
scale # cc d varies from 0.08" to 0.24" in steps of 
0.02", so the final Drizzled pixel sampling scales, 
9 Sl range from 0.04" to 0.12". Table 1 summarizes 
the four simulation sets . 

3.5. Weak Lensing Measurement 

We now turn to exploiting the simulated im- 
ages in order to recover the weak lensing signal 
and thereby determine the surface density of use- 
able galaxies under various assumed pixel sizes and 
exposure strategies. 

We locate galaxies in the noisy images using 
Source Extractor (Bertin & Arnouts 1996) con- 
figured to detect all objects as near to the noise 
threshold as possible via a Gaussian detection ker- 
nel matched to the known size of the chPSF at 
each pixel scale. The size of the detection kernels 
are therefore proportional to the chPSF FWHM 
values plotted in Figure 2. Source Extractor con- 
volves the images with this detection kernel, blur- 
ring features smaller than the size of the kernel, 
such as pixel-to-pixel photon counting noise, to 
avoid spurious detections. 

We adopted the shear measurement method by 
Rhodes, Refregier, & Groth (2000, hereafter RRG) 
to measure the shapes of detected galaxies. RRG 
was specifically developed with space-based weak 
lensing measurements in mind (Rhodes et al. 2000) 



and has undergone extensive tests on simulations 
(Leauthaud et al. 2007) and use on Hubble Space 
Telescope images (Rhodes et al. 2001; Refregier 
et al. 2002; Rhodes et al. 2004; Massey et al. 
2007a) to constrain cosmological parameters in- 
cluding er 8 , the normalization of the dark matter 
power spectrum. 

RRG is a modification of the KSB+ (Kaiser 
et al. 1995; Hoekstra et al. 1998) method which 
measures Gaussian- weighted multipole image mo- 
ments, 

Jii = J d 2 ew(9)I(e)0i0j, (12) 

where w is a Gaussian; i,j £ {x,y}, correspond- 
ing to the orthogonal image coordinates; and 
is chosen such that the weighted barycentcr is 
zero. RRG corrects the galaxy image moments 
for ePSF effects using the moments of the mea- 
sured ePSF. More advanced shape-measurement 
algorithms that are being developed, including the 
Shapelets-based method of Massey et al. (2007), 
instead model the chPSF. This requires higher res- 
olution data, but allows information known a pri- 
ori about a regular pixel response function to be 
included analytically, whether that be a fixed R sq 
or known pixel-to-pixel variations in the PRF. 

RRG forms a measure of ellipticity (in contrast 
to KSB+), only at the final stage, after correction 
of individual shape moments 



(Jx 



Jyy, 2J X y) 



J xx ~t~ Jyy 

A shear estimator is then formed, 



7 Q , 



(13) 



(14) 



where the shear susceptibility G is a scalar func- 
tion of higher order moments of the ensemble of 
galaxies. 

We then remove objects from the shear catalog 
in a similar manner to analyses of real data by 
Massey et al. (2007a); Leauthaud et al. (2007), 
and shown by the dashed lines in figure 5. We 
first eliminate the < 1% of galaxies that Source 
Extractor misclassified as stars, and also use 
the Source Extactor s/n outputs flux_auto and 
f luxerr_auto to cut galaxies with s/n < 10. 

We then calculate the 'RRG size', defined using 



7 



Table 1 
Summary of the simulations. 



Sim Set Perturbation a Plate Scale <7 ch b S C 



Sim 1 


Photosite size in /im 10"mm 1 


4^m = 0.04" 


#ccd 


Sim 1 Drizzled 






#ccd/2 


Sim 2 


Focal length (plate scale) 10"mm _1 x ^f^- 


4/zm = 0.04" x 


#ccd 


Sim 2 Drizzled 


^ccd/2 



a What is perturbed in order to change cc d- 

b Standard deviation of the Gaussian charge diffusion kernel as a function of perturbed CCD 
pixel scale. 

c The sample rate, as a function of perturbed CCD pixel scale, of the final images on which 
the weak lensing analyses are applied. 



the uncorrected quadrupole image moments as 



j-rrg 



'VV 



(15) 



We remove large galaxies with d IIg > 2000, 
which eliminates fewer than 1 detected object per 
arcmin 2 in our simulations. And via 



d IIg < 1.2y/<% + 0.4 = d c 



(16) 



we eliminate galaxies whose size is nearly equal 
to the measured ePSF. The 0.4 term also elimi- 
nates galaxies that are only a few pixels across, 
which have d rrg w 1. This term is negligible at 
small pixel scales, as Figure 3 shows, because the 
ePSF is large in pixel units. At large 9 cc d, where 
the ePSF size is nearly one pixel, the 0.4 term 
causes d cut — > 1 pixel, also seen in Figure 3. This 
cut ensures we measure shear on galaxies that are 
somewhat larger than both the PSF and the pixel 
size, which is the regime shape measurement is 
most reliable. The effect of the lower size cut on 
the data is shown with the dashed line below the 
ellipticity cut in Figure 5. We also test apply- 
ing a lower size limit of d cut = 1.2d* instead of 
Equation 16. This leaves n g (# C cd ^5 #o) unchanged 
within a few galaxies per arcmin 2 at most, but 
« g (#ccd > #o) increases by up to ~ 10/arcmin 2 . 
The result that smaller pixel scales always yield 
more galaxies, however, is robust. 



3.0 



2.5 



2.0 



1.0 



0.5 



0.0 

0.00 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

, 1.2 (dr e + 0.4) '/e. 
/ / in arcsec 




/ 

/ 

/ 

/ 

/ 

/ 

\ / / in arcsec 




\ \ / y ' 

\x 

- " "\. ~— _ 1.2 {d' s + 0.4) '/ B . 




in pixels 




— — _ 1.2 d- 




in pixels 







0.20 



0.15 



0.05 



0.00 



0.05 



0.10 0.15 
flood (arcsec) 



0.20 



Fig. 3. — Lower size limit c? cu t in pixels (solid, left 
axis) and arsec (dash, right axis) vs cc( j. (Note 
the left axis doesn't map horizontally to the right 
axis since # cc d is made to vary.) d* is the median 
d IIg measured on high s/n stars at each 8 cc d- This 
study applies cuts with the 0.4 term, which has 
miminal effect at small # C cd, but causes the size 
limit to tend toward about one pixel (dotted line) 
at large pixel scales. Plotted here are cuts made in 
the Sim 2 study, with size cuts in other simulation 
sets behaving similarly. 
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4. Results 



After discussing some properties of the simu- 
lated images, we explore as a function of pixel 
scale: the number of detected and "useful" galax- 
ies for weak lensing measurements n g ; the per- 
formance of shear recovery; the sample variance 
of shear estimators; and the predicted errors on 
the dark matter power spectrum. The number of 
galaxies for which shear can be measured is mainly 
a function of the angular size of galaxies relative 
to the size of the PSF. These results are there- 
fore likely to be independent of our adopted shear 
measurement method. However, the performance 
of shear recovery and the error on the power spec- 
trum will depend upon the sophistication of the 
chosen method, so those results are likely to im- 
prove before the launch of any future space-based 
mission. 

4.1. Image Properties 

As an illustration, Figure 4 shows the size and 
magnitude distributions of galaxies useful for lens- 
ing measurements in images from Sim 1 at the 
baseline pixel scale. The ordinate axis units are 
counts per arcmin 2 , and the data are measured on 
galaxies after we make cuts, as described in §4.2. 
Our images have a depth of about mag 27.7, but 
the faintest galaxies are not ultimately used. 

4.2. Surface Density of Sheared Galaxies 

The surface density of galaxies that survive the 
various cuts discussed above is plotted as the solid 
line in Figure 5. The uppermost dashed lines in 
each simulation set show the raw galaxy detections 
using Source Extractor 4 , averaged from all nine 
sheared images at each pixel scale in each Sim set. 
The surface density of galaxies that survive this 
cut are plotted in the next highest dashed line for 
each Sim set in Figure 5. 

As expected, the usable n g decreases monotoni- 
cally in all simulation sets, indicating that smaller 
# cc d always helps decrease counting noise in mean 
shear measurement for a fixed survey area. Stack- 
ing sub-pixel dithers fully recovers the surface den- 
sity of galaxies from Sim 1, despite larger the CCD 
pixel convolution at pixel scale 9 cc d — 29 s . Some 
galaxies are lost in the Sim 2 Drizzled data. Here 



4 See http://terapix.iap. fr/rubrique.php?id_rubrique=91/. 




0.01 I . I I . . I . . . I . . . I . I . I . . .1 III II I i 

20 22 24 26 28 0.1 1.0 

magnitude FWHM (arcsec) 



Fig. 4. — Surface density of galaxies as a function 
of magnitude (left panel) and size (right panel). 
These are the galaxies we use to measure the 
shear, which survive various cuts as described in 
§4.2, in Sim 1 images with # cc( j = 9$ (black) and 
#ccd = O.40o (dark gray), and in Sim 1 Drizzled 
images with # cc d = 29 (light gray). Magnitude 
bins are 0.25 wide, and size bins are 0.05 wide in 
log 10 space. Magnitude is taken from the Source 
Extractor mag_auto output, and FWHM from the 
f whm_image output. 
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Fig. 5. — The surface density of galaxies detected 
and usable for weak lcnsing, n g , along with counts 
after different cuts are applied to the data in turn, 
vs the pixel scale of the reduced images, S . The 
top panel plots data from Sim 1 (black lines) and 
Sim 1 Drizzled (gray lines) , and the bottom panel 
similarly for Sim 2 and Sim 2 Drizzled. Galaxies 
are detected with Source Extractor, and cuts are 
made based on measured ellipticity ("e cut"), size 
("d cut," see Equations 15 and 16), and s/n (final 
cut, giving "usable" galaxies). Error bars are cal- 
culated as the sample variance of n g between the 
different images at a given pixel scale. 



the constituent dithers, unlike the Sim 1 sets, have 
a larger underlying chPSF than the non-stacked 
counterparts. Taken together, these results indi- 
cate that n g worsens with increased charge diffu- 
sion, but is independent of the CCD pixel spacing 

#ccd- 

We estimate errors on the surface density of us- 
able galaxies, plotted with error bars in Figure 5, 
to be the sample variance of n g measurements be- 
tween each of the nine simulations per pixel scale. 
Errors are of order one galaxy per arcmin 2 for all 
data, which we consider to be negligible in our 
analysis. 

4.3. Accuracy of Weak Shear Recovery 

Figure 6 illustrates the accuracy with which we 
recover shear measurement with RRG for one typ- 
ical set of simulated images. 

We define 'accuracy' to be the closeness of 
the measured value to that originally input. In 
all cases, the shear recovery is well-fit by a lin- 
ear model thereby justifying the catalog cuts de- 
scribed in §3.5. Relaxing the cuts and using 
smaller or fainter objects introduces systematic ef- 
fects that cannot be corrected with RRG. Again in 
all cases, the y-intercepts ("additive shear resid- 
ual") of the shear recovery are consistent with 
zero and accordingly set explicitly so in subsequent 
analyses. 

On the other hand, the best- fit slopes ("mul- 
tiplicative shear bias") are systematically smaller 
than unity by roughly 20%. This effect has been 
known for some time as a limitation of KSB+ 
methods (eg Bacon ct al. 2001; Heymans et al. 
2006; Massey et al. 2007b), and it has been spec- 
ulated to arise as a result of pixelation. The pop- 
ulation bias has appeared to be robust to effects 
like galaxy morphological type, and we assume in 
all subsequent sections that it could, in practice, 
be determined to arbitrary accuracy using simula- 
tions. We therefore take 



bias 

b= + 1, 

7 



where, conventionally, 



bias = (7) — 7. 



(17) 



(18) 



to be the true value with negligible error, and cor- 
rect for the bias using values from the line fitted 
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DX = 0.10 arcsec 




-0.04 -0.02 0.00 0.02 0.04 

7i in 



Fig. 6. — Example of known, input shear vs 
measured shear, for one set of simulated images. 
Points are ( / y i ), measured from the data. Error 
bars are the standard error of shear estimators in 
each image. The dotted lines have slope 1, repre- 
senting perfect shear recovery, and the solid lines 
are linear least-squares fits to the data points. 
Also shown are the y-intercept a divided by the 
standard error of a, SE(a), and the slope b plus or 
minus the standard error of b, SE(b) — all outputs 
from the least-squares algorithm. Also shown are 
the x 2 values of the fit, which has 7 degrees of 
freedom. 



as a function of pixel scale in figure 7. We call b 
the bias for brevity. 

4.4. Precision of Weak Shear Recovery 

Weak lensing induces changes of only a few per- 
cent in galaxy ellipticities, but the RMS of the 
cllipticity distribution of faint galaxy populations 
are about 30%, (e.g. Leauthaud et al. 2007). Fig- 
ure 8 shows the standard deviation of ellipticities 
cr e and shear estimators ct 7 , where a e is defined 
analogously to 

CT 7=<+<' ( 19 ) 

where 

<H(7,-<7«» 2 >. (20) 

The measured RMS cllipticity increases mono- 
tonically with cc d, meaning larger pixels always 
make shape measurements noisier, an undesirable 
effect for weak lensing analysis. The ePSF "di- 
lutes" the shear signal in galaxies, and RRG is 
shown here to reverse this dilution at all pixel 
scales by decreasing the RMS of shapes, as ex- 
pected. The near coincidence of the non-Drizzled 
and Drizzled <r 7 lines suggests that the RMS shear 
depends most on the CCD pixel spacing 9 CC< \ 1 and 
whether or not ideal dcintcrlacing is performed. 
This is, in fact, the ideal situation: the best im- 
age resolution we can have is fundamentally lim- 
ited by (in addition to the diffraction) the CCD 
pixel response — that is, the ePSF. Ideal deinter- 
lacing increases the sample rate of a given object, 
which is why a e decreases when Drizzling. This 
plot shows that RRG recovers the underlying, ch- 
PSF convolved shape information from the ensem- 
ble even without Drizzling. 

By assuming uncorrelated shapes in galaxy 
populations (for a discussion, see Hirata & Sel- 
jak 2004), this noise can be reduced by measuring 
many galaxies and applying Poisson statistics. An 
ability to use more galaxies per arcmin 2 , n g , im- 
proves the precision of shear recovery, which is 
best quantified by the sample variance of shear es- 
timators, or its square root a^/y/n. This is shown 
as a function of pixel scale in Figure 9. This re- 
iterates our result that smaller pixel scales always 
improve the precision of weak shear estimation. 
However, smaller pixel scale face diminishing re- 
turns: the loss of precision when increasing the 
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Fig. 7. — The RRG shear estimator bias b, and 
deviation from the least-squares line fitted to each 
component Ab, as a function of pixel scale 6> cc d- 
Squars indicate the 71 component of shear, rela- 
tive to the pixel grid, and diamonds the 72 compo- 
nent. Error bars show the standard error on slope 
measurement obtained from the least-squares fits 
to data like those show in figure 6. The bias in 
each component of shear is different, and varies 
with pixel scales. Extrapolating the pixel scale to 
zero appears to give a bias of about 0.8, so pixe- 
lation is likely not the cause of the bias. Drizzling 
to smaller pixel scales does not remove the bias or 
the trend toward 0.8. 
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Fig. 8. — RMS galaxy shapes in single expo- 
sure (black lines) and Drizzled (gray lines) sim- 
ulations. T]he_ilMS of the shear is calculated as 
cr 7 = \j<^\ + , where a 7 i is the RMS of the 
bias-corrected shear estimators. Errors are calcu- 
lated as the sample variance of the RMS shapes 
between different images at a given pixel scale. 
The RMS of the galaxy ellipticities are calculated 
similarly from uncorrected ellipticities. 
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pixel scale by some amount is greater in mag- 
nitude than the gain in precision of decreasing 
the pixel scale by the same amount. Drizzling 
decreases the error somewhat, thanks entirely to 
the behavior of n g with cc( j. The shear errors 
from the Sim 2 sets degrade faster with increased 
9 cc d- This indicates that charge diffusion, which 
dominates larger 9 cc d in the Sim 2 sets, should be 
minimized in addition to # C cd- 

One potential concern is the observation by 
Kaiser (2000) that the distribution of shear esti- 
mators in practice is not Gaussian, and that its 
extended wings may even make its second moment 
infinite. This is not surprising because the shear 
estimator involves a ratio of two noisy quantities. 
Figure 10 shows the distribution of RRG shear es- 
timators measured in all 9 images at different pixel 
scales, and Figure 11 shows the kurtosis of those 
distributions. The distributions are indeed non- 
Gaussian. However, there does not appear to be 
any significant degradation in this sense at larger 
pixel scales. The skewness of the distributions is 
also consistent with zero at all pixel scales. 

4.5. Cosmological Implications 

As we have already seen, the quality of shear 
measurements is always improved with small pixel 
scales. However, for a mission with a fixed lifetime, 
larger pixels would allow a linear increase in the 
total survey volume, and a corresponding decrease 
in sample (or "cosmic") variance errors. These 
two effects combine in a measurement of the dark 
matter power spectrum from cosmic shear, which 
would have a total statistical error (c/ Refregier 
ct al. 2004) 



AC e = 




(21) 



where d is the power, and t is a multipole. The 
fraction of sky surveyed is 



/sky {0, 



ccd ) 



r P ccd 



(22) 



where /o is the baseline fraction, 9q the baseline 
pixel scale, and / s k y at the largest simulated pixel 
scale is assumed to be still smaller than the observ- 
able sky. This introduces a tension between weak 
lensing precision and cosmic variance. We use the 
overall error AC; as our final figure of merit. 
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Fig. 9. — The standard error of the shear, cr n 
in the single exposure (black lines) and Drizzled 
(gray lines) simulations. Error bars are calculated 
as the sample variance of the standard error be- 
tween the different images at a given pixel scale. 
This quantity combines the method-independent 
n g and the more method-dependent er 7 into an 
estimate of the precision of our shear estimation. 
This is the figure of merit showing how precisely 
we can measure shear as a function of flccd, if the 
survey area is fixed. 
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Fig. 10. — The distribution of shear estimators 
71 obtained from Sim 1 images with # cc d = 0.10" 
(black) and cc( j = 0.04" (dark gray), and Sim 
1 Drizzled images with # cc d = 0.20" (light gray) 
pixel scales. Each sample comes from all images 
at the quoted pixel scales, which have different in- 
put shears applied. The mean of the input shears 
at each pixel scale, however, is zero, so these dis- 
tributions are themselves centered around zero. 



Fig. 11. — Kurtosis of shear estimators minus 3 
MS (X function of pixel scale. The kurtosis is calcu- 
lated as ((7» — (7i)) 4 ) / (Tj. , where er 7i is the variance 
of 74. The 7^ are measured from all 9 images at 
each # C cd- 3 is the kurtosis of a Gaussian distri- 
bution. All data here are greater than 0, showing 
that the distribution of measured shear estimators 
has a sharper peak and more extended wings than 
a Gaussian. 



14 




5 ........................ , 
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flood (arcseo) « / sk y 1/2 

Fig. 12. — The contribution of shear sample vari- 
ance to the predicted error on the matter power 
spectrum in single exposure (black lines) and Driz- 
zled (gray lines) simulations. This assumes the 
number of pixels in the focal plane is fixed when 
perturbing 9 cc d, such that the survey area scales 
as in Equation (22). The ordinate axis scale is 
log 10 , normalized to the Sim 1 error value at the 
baseline pixel scale, 9 cc d = 0.10". 



The measured values of the second term (i.e. 
in the limit of small Ci) are shown in Figure 12. 
The ordinate axis uses a logarithmic scale, in units 
of the baseline survey because the absolute val- 
ues are heavily dependent on the the survey area, 
multipole £, number of galaxies, etc. Cosmic vari- 
ance dominates errors at small # cc d, and increased 
sample variance in shear estimators takes effect at 
large # cc d- The error flattens out at larger pixel 
scales, affording some freedom in choose the pixel 
scale, in which case the smaller pixel scales are 
clearly preferred for reasons of caution. 

Our idealized Drizzling decreases the error, 
thanks entirely to the behavior of n g with C cd> 
indicating larger CCD pixel scales are acceptable 
if perfect deinterlacing is possible and charge dif- 
fusion remains fixed. It is important to note that 
larger # cc d degrades n g , RMS shapes, and ulti- 
mately shear errors. As expected, charge diffusion 
dominates large at # cc d in Sim 2 sets. 

Increasing the survey lifetime while fixing the 
depth would increase the nominal survey area /o, 
and thus shifts all curves down uniformly on a 
log 10 scale; however, this may increase the range 
of angular scales probed (£) and thus change the 
overall behavior of Equation (21) with 9 cc d- 

5. Discussion 

To explore the effects of pixclation on galaxy 
shape measurement, we have realistically simu- 
lated weakly sheared galaxies at a range of pixel 
scales. We have examined two different ways in 
which a future weak lensing space mission could 
alter the pixel scale, and both non-dithered and 
dithered exposure strategies, for a total of four 
different simulation sets. We have then emulated 
the process of shear measurement that is applied 
to real data. We have finally combined the ob- 
served surface density of galaxies, statistics of the 
shear estimators, and our privileged knowledge of 
the true input signal, to arrive at figures of merit 
for the pixel scale. 

We find that, as expected, smaller pixel scales 
consistently give improve the quality of shear re- 
covery and tighten constraints on the dark matter 
power spectrum for a hypothetical survey of fixed 
area. Ideal deinterlacing gives a further improve- 
ment. On the other hand, when the survey area is 
permitted to change with the pixel scale according 
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to larger or smaller detectors, we find that larger 
pixel scales minimize statistical errors on the mea- 
sured matter power spectrum because the survey 
area increases as 0^ cd - These errors flatten out 
somewhat above 6> cc d = 0.09". Considering both 
situations of fixed and variable survey area, the 
best # cc d would be the smallest allowable by the 
projected dark matter constraints, which we find 
to be about 0.75-0.80 the FWHM of the charge 
diffusion convolved diffraction pattern, chPSF. 

We used a current-generation shear measure- 
ment method as the basis of our study — a snap- 
shot of available technology. Better methods are 
certainly needed to fully exploit the ambitious fu- 
ture surveys now being planned. We can specu- 
late that these may either be better able to cope 
with poor resolution, or (more likely, since infor- 
mation is irrevocably lost during pixelation) re- 
quire smaller pixels to overcome systematic floors 
revealed by the lowering of statistical errors. This 
is being suggested by the variable bias in current 
results, which is now well documented but poorly 
understood. It even has to be argued whether the 
best use of an expensive space mission would be to 
minimize statistical errors on an isolated measure- 
ment of the matter power spectrum. By imaging a 
smaller region, but with a higher density of useable 
galaxies, and smaller errors, a mission could alter- 
natively be used to obtain the phase information 
needed for maps or for higher order correlation 
functions; or even to calibrate the shape measure- 
ment of larger, ground-based surveys. As shown 
by our results, this approach would prefer smaller 
pixel scales. 

Two important simplifications were imposed on 
our pipeline. Firstly, we allowed no temporal or 
spatial variation in the PSF, and a created a com- 
fortably large number of fake stars to characterize 
the PSF shape. As demonstrated in Rhodes et al. 
(2007), pixelation especially adds noise to peaky 
objects like a diffraction-limited PSF. A typical 
survey region is likely to lie at high galactic lati- 
tude. If sufficiently bright stars cannot be imaged 
within the time taken for the PSF to vary, noise in 
the measurement of PSF shapes (which we have 
not considered) could potentially dominate that 
in galaxy shapes. Secondly, real dither strategies 
never provide perfect interlacing. For example, op- 
tical distortions differentially alter the spacing of 
the pixel grid in different places, effectively caus- 



ing 9 cc d to be a function of position in the focal 
plane. Consequently, Drizzling real images cor- 
relates noise between adjacent pixels because it 
must average nearby pixel values. Correlated noise 
hinders both object detection and shape measure- 
ment, so Drizzling is detrimental to weak lensing. 
The dithering implemented here is therefore ide- 
alized, and provides a best-case scenario. 

With these caveats in mind, our approach has 
provided a practical analysis that is achievable 
with existing methodology. It is complementary to 
studies starting from analytic first principles and 
assuming the existence of a perfect shape mea- 
surement method. In practice, our result on the 
optimum pixel scale should sensibly be considered 
as an upper limit, pending future developments in 
shape measurement methodology. 
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